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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Sysiems  Command  (AFSC).  The  results  were  obtained  by 
Calspan  Field  Services,  Inc.,  operating  contractor  for  the  Aerospace  Flight  Dynamics  testing 
crforl  at  the  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee.  The  Air  Force  project 
manager  was  Mr.  Marshall  Kingery.  This  work  was  done  under  AEDC  project  numbers 
V32S-P4  and  V-45G.  The  manuscript  was  submitted  for  publication  on  June  10,  1981. 
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1.0  INTRODUCTION 

An  appreciable  capability  exists  in  the  use  of  curve-fitting  computer  programs  for 
obtaining  static  stability  constants  from  motion  histories  of  bodies  having  forces  and 
moments  that  are  nonlinear  with  angle  of  attack.  Such  motion  histories  can  be  obtained  in 
either  aeroballistic  ranges  or  wind  tunnels.  However,  as  demonstrated  in  Ref.  1,  there  are 
situations  for  which  it  is  advantageous  to  obtain  fits  of  “nonlinear  motion  histories”*  by 
using  curve-fitting  programs  in  which  only  linear  terms  of  the  motion  equations  are  used, 
thereby  providing  linear  fits  of  the  motion  histories.  When  only  linear  terms  are  used  in 
fitting  nonlinear  motion  histories,  mean  static  moment  and  force  parameters  are  obtained. 
These  are  referred  to  as  “effective”  stability  derivatives  which  are  designated  and  CNa  , 
respectively,  in  this  report.  For  simplicity,  only  the  effective  moment  derivative  will  be 
discussed  here,  but  comments  made  will  also  apply  to  the  corresponding  effective  force 
derivative. 

Because  the  effective  moment  derivative,  C;^ ,  is  measured  over  an  interval  of  a  motion 
history  and  hence  corresponds  to  a  band  of  motion  amplitude,  defining  the  relationship  of 
Cm0  with  the  corresponding  nonlinear  moment  coefficient,  Cmt,  has  been  a  problem  area  in 
the  past.  Cm[  is  the  desired  parameter  from  experimental  tests  of  nonlinear  configurations, 
and  measurements  of  Cmi  are  usable  when  they  are  defined  either  as  a  function  of  the  angle 
of  attack  or  tabulated  versus  corresponding  angle-of-attack  values. 

An  approximate  procedure  for  relating  to  the  corresponding  nonlinear  moment 
coefficient,  using  an  effective  angle-of-aftack  concept,  was  presented  recently  in  Ref.  1 .  This 
procedure  was  shown  to  be  very  useful  in  the  analysis  of  the  nonlinear  experimental  data  of 
Ref.  1.  The  purpose  of  this  report  is  to  discuss,  with  the  use  of  the  effective  angle-of-attack 
concept  of  Ref.  1 ,  the  use  of  effective  moment  derivative  measurements  obtained  by  various 
investigators  at  different  lest  conditions,  (e.g.  Refs.  2  to  8).  Other  approximate  techniques 
such  as  those  described  in  Refs.  9  and  10  have  been  published  that  concern  relating  a 
measured  effective  moment  derivative  to  its  corresponding  nonlinear  moment  coefficient. 
Those  techniques  involved  more  general  approaches  to  the  problem  and  do  not  permit  an 
assessment  of  the  restrictions  related  to  the  direct  utilization  of  experimental  effective 
derivatives. 


‘Motion  histories  corresponding  to  bodies  having  nonlinear  force  and  moment  characteristics. 
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2.0  DISCUSSION 
2.1  EFFECTIVE  ANGLE-OF-ATTACK  CONCEPT 

The  procedure  presented  in  Ref.  1  for  approximating  the  relationship  of  an  effective 
moment  derivative,  Cjt  T  to  its  corresponding  nonlinear  moment  coefficient  is  briefly 
described  in  the  following  paragraphs. 

For  an  ideal  linear  configuration,  the  moment  coefficient  can  be  written 

C  “  C  •  a  rn 

mt  ma  (U 

where  Cm  a  is  a  constant  for  different  values  of  a.  If  the  oscillatory  motion  histories  of  such  a 
linear  configuration  are  fitted  with  a  computer  program  (using  only  the  linear  terms  of  the 
program),  Cllla  of  Eq.  (I)  is  one  of  the  parameters  defined.  The  amplitude  of  the  motion 
interval  fitted  is  of  no  concern  because  Cm<i  is  a  constant  and  Cm,  can  be  defined  for  the  a 
values  of  interest. 


For  a  configuration  having  a  nonlinear  moment  that  can  be  described  with  a  cubic 
nonlinearity  (a  cubic  nonlinearity  is  normally  assumed  in  the  fitting  technique),  then 


C_  ’C. 


a  +  C.. 


't  a 


■3  a 


(2) 


Here,  a  is  used  as  the  angle  of  attack  (small  angles)  or  sine  of  the  angle  of  attack.  If  a  motion 
history  of  such  a  configuration  is  fitted  as  above,  using  only  the  linear  terms,  then  the 
motion-fitting  program  will  define  a  quasi-linear  moment  parameter,  Cj^a.  The  usefulness 
of  C„u  values,  as  noted  in  Section  1.0,  has  been  limited  in  the  past  by  a  lack  of  techniques 
for  utilization  of  corresponding  amplitude  information  available  from  the  motion  fit  to 
deduce  nonlinear  moment  characteristics.  For  example,  if  C^q  measurements  were  obtained 
from  motion  fits  of  a  given  configuration  at  different  amplitude  levels  and  those 
measurements  were  approximately  constant,  then  such  measurements  would  have  confirmed 
that  the  moment  variation  for  the  configuration  was  approximately  linear.  However,  if 
different  values  were  obtained  for  Cma,  then  one  would  know  only  that  the  configuration 
had  a  nonlinear  moment  and  that  Cm[,  which  is  of  primary  concern  in  a  nonlinear  system, 
could  not  be  defined  directly. 


If  a  value  has  been  determined  from  a  fit  of  an  interval  of  motion  history  of  a 
vehicle  that,  in  fact,  has  a  nonlinear  moment  coefficient  defined  by  Eq.  (2),  then  it  follows 
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that  there  is  some  effective  angle  of  attack,  ae,  of  the  motion  interval  fitted  that  corresponds 
to  the  such  that  Cm(  of  Eq.  (2)  can  be  written 


(3) 


The  obvious  problem  here  is  in  relating  ae  to  the  amplitude  of  the  motion  interval  fitted.  To 
accomplish  this,  the  ratio  (ac/0m)  is  used.  Here,  0m  is  the  average  of  (a)  the  midrange 
amplitude  of  the  fitted  interval  of  motion,  and  (b)  the  mean  of  the  initial  and  final  motion 
amplitudes  of  the  fitted  interval.  Further,  0m  corresponds  to  the  major  axis  in  the  case  of  an 
elliptic  angle-of-attack  motion  pattern.  If  ac  can  be  defined  adequately,  then  one  point  on 
the  nonlinear  Cmt  curve  can  be  defined  for  each  fitted  Csra  value  without  the  problem  of 
evaluating  both  Cm)a  and  Cm3(>  of  Eq.  (2). 


The  results  of  a  study  of  the  ae  concept  presented  in  Ref.  I  for  realistic  free-flight 
conditions  indicate  that  the  (crc/0m)  ratio  tends  to  be  insensitive  to  the  magnitude  of  the 
nonlinearity,  to  the  type  of  motion  pattern,  and  to  the  amplitude  level  involved  for  a  fairly 
wide  range  of  test  conditions.  The  study  utilized  a  series  of  computer-generated  nonlinear 
motion  histories  for  given  Cmin  and  Cm3c(  values  which  were  fitted  (linear  fits  utilizing  the 
correct  roll  rates)  to  provide  corresponding  Cma  values.  The  ae/0m  ratio  was  shown  to  be 
nearly  constant  with  a  mean  value  of  0.87,  and  the  corresponding  Cmi  values  can  be  obtained 
quite  well  using  the  ac  concept  from  nonlinear  free-flight  motion  histories,  given  the 
following  restrictions  that  (1)  the  test  configuration  is  statically  stable  at  zero  angle  of 
attack,  (2)  the  magnitude  of  the  nonlinear  Cm(  curve  corresponding  to  the  maximum 
amplitude  of  the  motion  interval  being  fitted  is  not  less  than  about  0.9  of  the  Cm|  peak  value 
of  the  interval,  (3)  the  ellipticity*  of  the  motion  pattern  is  less  than  about  0.5,  and  (4)  the 
model  roll  rate  is  small.  The  relationship  between  C^a,  Cm|,  ae,  and  0m  is  shown  below. 


*Raiio  of  the  lengths  of  the  minor  axis  to  the  major  axis  of  the  elliptic  motion  pattern  experienced 
by  a  body. 


7 


A  E  DC-TR-81  -1  2 


2.2  LARGE  AMPLITUDE  MODE  OF  TESTING 

To  demonstrate  the  limited  use  that  has  been  made  of  large  amplitude  moment  data,  the 
measurements  of  of  Fig.  2  of  Ref.  2  have  been  replotted  in  Fig.  la.  These  are 
measurements  for  cones  with  different  nose  bluntness  ratios,  as  obtained  with  a  large 
amplitude  gas  bearing  pivot  system  (planar  motion);  they  are  representative  of 
measurements  of  a  series  of  large  amplitude  tests  that  have  been  conducted.  Corresponding 
predicted  moment  coefficient  curves  (Ref.  1 1)  for  the  two  blunt  cones  are  shown  in  Fig.  lb. 
These  curves  have  appreciable  nonlinearities. 


svj. 

O  0.001  O.POO  0.600  0.0009  2.-12  )  In  Ref.  2 

□  0.151  0.542  0,600  0.0015  3.19  L  -  Virtual  Cone  Length 

O  0.300  0.465  0.600  O.OOL2  2.72  !  i.  -  Actual  Cone  Length 

M  =  10 

SO 

-  Unsteady  Flow  Field  Theory 

Note:  of  Ref.  2  is  same  as  C—  of  present  report. 


a.  Large  amplitude  measurements  (Fig.  2,  Ref.  2) 
Figure  1.  Moment  daita  for  10-deg  semiangle  cones. 
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As  previously  discussed,  in  testing  nonlinear  configurations  it  is  the  Cm|  values  that  are  of 
primary  concern.  Reference  2  notes  that  the  measurements  were  not  constant  with 
increasing  amplitude  for  the  two  blunt  cones,  indicating  that  the  corresponding  moment 
coefficient  curves  would  not  be  linear.  This  was  an  appropriate  comment  because  a  direct 
relationship  between  a  measurement  and  its  corresponding  Cm(  value  was  not  available 
at  that  time. 


b.  Predicted  curves.  Ref.  11  (M  =  10) 
Figure  1.  Concluded. 


In  Fig.  2  the  predicted  moment  curves  of  Fig.  lb  are  shown  again,  along  with 
experimental  Cm|  values  computed  using  C„a  measurements  of  Fig.  la  at  amplitudes  of  4 
and  8  deg  and  utilizing  the  effective  angle-of-attack  concept  [with  («e/0m)  =  0.87].  The 
agreement  of  the  predicted  curves  with  the  experimental  Cffl[  values  in  Fig.  2  provides 
additional  confidence  in  the  use  of  the  effective  angle-of-attack  concept  for  obtaining 
nonlinear  moment  coefficients  from  measurements.  A  well-defined  experimental  Cm[ 
curve  can  be  obtained  by  using  additional  measurements  from  Fig.  la.  The  Cma 
measurements  of  Ref.  2  would  be  expected  to  be  of  high  quality  because  they  were  obtained 
with  an  air  bearing  pivot  system  which  would  have  had  a  trivially  small  tare  stiffness 
parameter.  Normally,  the  measuring  error  in  measurements  in  carefully  designed  tests 
can  be  kept  to  within  one  percent. 
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Figure  2.  Comparison  of  Cmt  values  computed  from  Csa 
data  of  Fig.  la  with  predicted  curves  of  Fig.  1b. 

2.3  SMALL  AMPLITUDE  MODE  OF  TESTING 


2.3.1  Fora,  =  0 

In  this  mode  of  testing,  previous  investigators  (see,  for  example,  Refs.  3  and  4)  have 
made  C^a  measurements  at  amplitude  levels  between  ±  1  and  +3  deg.  Because  of  the  lack 
of  information  in  the  use  of  measurements,  it  has  usually  been  assumed  that  (C^  ) 

is  equal  to  (3Cm,/3a)a=0.  The  validity  of  such  an  assumption  is  questionable  as 
demonstrated  by  the  differences  between  (3C 0  and  (C^)  W|  =  0  values  shown  in  Fig. 
3a.  The  curve  of  Fig.  3a  corresponds  to  a  configuration  having  a  representative  nonlinear 
moment  curve  (cubic  nonlinearity  with  Cm]o  =  -0.542  and  Cmja  =  +  16.5,  as  shown  in 
the  sketch  on  p.  11).  This  same  assumed  moment  curve  is  used  in  assessing  measurements 
made  under  other  test  conditions,  as  discussed  in  later  sections  of  this  report. 
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Angle  of  .Attack,  deg 


The  curve  in  Fig.  3a  was  obtained  from  a  computer  study  in  which  linear  fits  of  generated 
nonlinear  motion  histories  were  obtained.  Figure  3a  shows  that  the  difference  between 
(Cm„)ct  =  0(dClt)/3a)a  =  0  increases  with  the  amplitude  level  of  the  lest  and  is  about  seven 
percent  at  an  amplitude  of  ±3  deg.  Hence,  to  maximize  the  usefulness  of  C™  measure- 

111  QT 

ments,  it  is  important  to  always  minimize  the  amplitude  level  used  in  experimental  tests 
involving  nonlinear  configuration.  This  comment  follows  in  consideration  that  testing  can 
involve  configurations  having  nonlinearities  appreciably  greater  than  that  of  the 
demonstration  Cni|  curve  shown  above.  Any  nontrivial  measuring  error  associated  with  the 
particular  test  equipment  used  could  be  an  additional  factor  in  selecting  the  optimum 
minimum  test  amplitude. 
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Mote:  /C-  \  values  obtained  from  linear  fits  of 

\  a/at=0 

nonlinear  motion  histories  —  cubic  function 
with  C  =  -  0.542  and  C  =  +  16.5 
mla  m3« 
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a.  at  =  0 

Figure  3.  Comparisons  of  Cs;a  and  3Cmi/9o  values. 


2.3.2  For  cr,#0 

In  this  mode  of  testing,  investigators  in  the  past  (see,  for  example,  Refs.  4  through  8) 
have  used  a  measured  as  being  equal  to  (<3Cm /3a)a=U[  again  because  of  a  lack  of 

amplitude  information  for  corresponding  measurements.  Differences  between 

ancl  (®Cm/^a!)o[  =  al  that  can  expected  for  a  representative  nonlinearity 
(same  as  used  above)  are  shown  in  Fig.  3b.  The  curves  in  Fig.  3b  were  determined  again  in  a 
computer  study  in  which  linear  fits  of  generated  nonlinear  motion  histories  were  obtained. 
In  this  case  a  pivot  system  with  a  tare  restoring  moment  (having  linear  characteristics)  was 
simulated  using  (Cmu)pivm  =  -3.0;  this  permitted  providing  a  trim  moment  to  the  model 
and,  in  turn,  a  trim  angle  when  the  support  sting  was  set  at  a  nonzero  angle  of  attack.  The 
curves  of  Fig.  3b  indicate  that  the  difference  between  (C^  ^0and  (dCm  /da)a^a  in  small 
amplitude  testing  increases  appreciably  as  the  trim  angle  is  increased.  Again  in  such  tests  it  is 
important  to  minimize  the  amplitude  level  used.  It  is  significant  to  observe  for  a  body  having 
rotational  symmetry  that  a  nonlinear  Cmt  curve  corresponding  to  small  amplitude  motion 
about  a,*0  is  asymmetric,  and  the  use  of  (ctc/6m)  =  0.87,  derived  for  oscillatory  motion 
about  at  =  0,  is  not  valid.  Further,  one  can  expect  (0^/0,^  10  an  appreciably  more 

sensitive  parameter  than  («c/0m)u  Because  of  the  differences  between  (Cma)„  and 
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(dCm/da)a^,  which  can  be  significant,  one  must  exercise  care  in  the  use  of  a  procedure 
utilized  in  the  past  in  which  a  Cm|  curve  is  determined  from  an  integration  involving  the 
corresponding  (C^  )a  measurements. 


Amplitude  of  Motion,  deft 


b.  at  =£  0 

Figure  3.  Concluded. 


2.4  OUT-OF-PLANE  MEASUREMENTS 


In  this  mode  of  testing  (see,  for  example,  Refs.  6  through  8)  the  model  has  a  trim  angle  of 
attack  in  one  plane  and  small  amplitude  motion  in  the  opposite  plane.  As  an  example,  a 
typical  set  of  test  conditions  could  be  a,  *  0  and  0  =  ±2  deg.  For  this  mode  of  testing,  it  is 
useful  to  examine  the  corresponding  moment  equations  involved  for  a  body  having 
rotational  symmetry.  For  such  a  body  the  equation  for  the  total  moment  coefficient,  Cm(, 
corresponding  to  a  cubic  nonlinearity  was  defined  previously  (Eq.  2): 
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Components  of  a  are  defined  below. 


Note:  y  and  2:  are  two  axes  of  an  orthogonal  axis  system 

that  has  its  x-axis  aligned  with  the  longitudinal 
axis  of  the  body. 


Here  the  orientation  of  the  total  angle-of-attack  plane  relative  to  the  orthogonal  stability 
axes  (y  and  z)  is  a  function  of  the  angle  a,  which  is  defined  by  the  velocity  components  v  and 
w  such  that  tan  a  =  v/w.  For  stability,  the  total  moment  coefficient  Cmt,  is  oriented  in  the 
sketch  consistent  with  the  given  direction  of  a.  Note  that  in  conventional  motion  equations, 
Cm  and  Cn  are  the  desired  components  of  Cmt  relative  to  the  y  and  z  axes  and  will  have 
opposite  signs.  They  can  be  written 


C  =  C  •  sin  <7 

n  ri  . 


(4) 


C 


m 


C  •  cos  a 
m. 


(5) 


From  symmetry  considerations  for  a  body  with  rotational  symmetry  where  Cm(  is 
independent  of  a  and  considering  the  sign  difference  between  Cm  and  Cn  noted  above,  the 
equation  constants  of  Eq.  (1)  can  be  written  in  terms  of  the  conventional  stability 
coefficients  as 


'la 


=  C 


Dij8 


-  -G_ 


la 


(6) 
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and 


C;  =  c  =  -C 

J,n3a  n3^3  m3tt 

Equations  (4)  and  (5)  expressed  as  functions  of  the  velocity  ratios  can  be  written 


(7) 


C 


n 


(8) 


=  C. 


n 


1/8 


=  C„  Cn 

“IjQ  n3/3  n3/S 


and 


C 


m 


=  c 


m 


+  c. 


(9) 


+  C„ 

11 3<i 


The  partial  derivaties  of  Cn  and  Cra  [Eqs.  (8)  and  (9)]  relative  to  /3  and  a,  respectively,  are 


+  3C 

n3/8 


(10) 


and 


era 


=  c. 


'la 


+  3C. 


1:1 35 


-2  r  o2 

a-  +  C  j8 

m3a 


(11) 


These  equations  are  instructive  in  that  Eqs.  (8)  and  (10)  demonstrate  that,  even  for  a  body 
with  rotational  symmetry,  both  the  moment  coefficient,  Cn,  and  its  corresponding  partial 
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derivative,  (3Cn/3j3),  are  dependent  on  the  angular  motion  of  the  model  in  the  opposite 
plane.  A  similar  comment  could  be  made  regarding  Eqs.  (9)  and  (11). 

It  follows  from  Eqs.  (10)  and  (1 1)  that  the  in-plane  partial  derivative  (5Cm/3a)(/3=oi-=-i) 
can  differ  appreciably  from  the  out-of-plane  partial  derivative  (3Cn/a/3),3=0i-  =  -t), 
and  an  example  of  such  possible  differences  is  shown  in  Fig.  4.  The  effective 
stiffness  parameter,  (Cii0)(jS,-oIn-ai)i  measured  in  out-of-plane  tests  can  differ  from 
(3Cn/dd)((a =0,^=0,)  and  differences  that  can  be  expected  for  the  same  assumed  representative 
nonlinearity,  utilized  previously,  are  shown  in  Fig.  5.  The  curves  in  Fig.  5  again  demonstrate 
the  importance  of  always  minimizing  the  motion  amplitude  that  is  used  in  such  tests. 

From  these  comments  concerning  out-of-plane  measurements,  it  follows  that  an 
investigator  can  have  a  dilemma  in  assessing  the  usefulness  of  (Cn#J)(j3|  =0>-  =  -t)  type 
measurements  in  the  general  case  of  an  asymmetric  body. 


Note:  Curves  computed  for  cubic  moment  with 
C  =  -  0.542,  C  =  +  16.5 


Figure  4.  In-plane  and  out-of-plane  moment  data. 
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Figure  5.  Comparisons  of  and  0C„/3/3  values. 


CONCLUDING  REMARKS 

A  study  was  made  of  the  utilization  of  measurements  of  the  effective  moment  derivative 
as  obtained  in  small  amplitude,  large  amplitude,  and  out-of-plane  modes  of  testing.  This 
study  made  use  of  the  approximate  relationship  between  a  measurement  and  its  parent 
nonlinear  moment  coefficient,  Cm[  ,  as  defined  in  Ref.  1.  Results  of  the  study  indicate  the 
following: 

1.  A  nonlinear  moment  curve  can  be  defined  very  well  with  Cm„  measurements  and 
corresponding  effective  angles  of  attack  obtained  in  large  amplitude  tests. 

2.  In  small  amplitude  testing,  it  is  important  always  to  minimize  the  motion  amplitude 
if  the  measured  effective  parameter,  C-q  ,  is  to  adequately  approximate  the  desired 
corresponding  dC^/da  parameter.  This  is  true  whether  the  small  amplitude  tests  are 
for  a  zero  or  non-zero  trim  angle  of  attack.  A  nontrivial  measuring  error  associated 
with  the  particular  test  equipment  used  could  be  an  additional  factor  in  selecting  the 
optimum  minimum  test  amplitude. 

3.  One  must  exercise  care  in  assessing  the  usefulness  of  out-of-plane  measurements 
in  the  general  case  of  an  asymmetric  body. 
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NOMENCLATURE 


C 


mt 


c 


m 


c 


n 


Cm,, 


Moment  coefficient  corresponding  to  the  total  angle-of-attack  plane 

Moment  coefficient  about  the  y-axis  (see  page  14) 

Moment  coefficient  about  the  z-axis  (see  page  14) 

Moment  stability  derivative 

Effective  moment  parameter  obtained  in  a  linear  fit  of  a  nonlinear  moment 
history 


C  C 

C  -  C  - 

c  c 

a 


V 

v 

w 

a,  0 


Coefficients  for  a  cubic  moment  variation  associated  with  the  total  angle-of- 
attack  plane  [see  Eq.  (2)] 

Coefficients  for  a  cubic  moment  variation  —  associated  with  y  and  z 
components  of  the  angular  motion  [see  Eqs.  (7)  and  (8)] 

Total  angle-of-attack  (small  angles)  or  sine  of  the  total  angle  of  attack  (large 
angles) 

Total  velocity  vector 

Velocity  component  along  the  y-axis 

Velocity  component  along  the  z-axis 

Components  of  a  ' 


Effective  angle  of  attack,  [see  Eq.  (3)] 
at,  at,/3,  Trim  angles 

0m  Average  of  (a)  the  midrange  amplitude  of  the  fitted  interval  of  motion  and 

(b)  the  mean  of  the  initial  and  final  motion  amplitudes  of  the  fitted  interval 
Note:  All  moment  data  are  computed  using  the  model  diameter  as  the  reference  length  and 
the  model  center  of  gravity  as  the  moment  reference. 
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